Ground State Phase Diagrams of Frustrated Spin-S' XX Z Chains: 

Chiral Ordered Phases 
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The ground-state phase diagram of the frustrated spin-/? XX Z chain with the competing nearest- 
and next-nearest-neighbor antiferromagnetic couplings is studied numerically by using the density- 
matrix renormalization-group method for the cases of S = 1/2, 3/2, and 2. We are paricularly 
interested in the possible gapless and gapped chiral phases, in which the chirality k; = SfSf +1 — 
Sf Sf +1 exhibits a finite long-range order whereas the spin correlation decays either algebraically or 
exponentially. We show that the gapless chiral phase appears in a broad region of the phase diagram 
for general S. By contrast, the gapped chiral phase is found for integer S in a narrow region of the 
phase diagram, while it has not been identified for half-odd integer S within our numerical accuracy. 
By combining the results with our previous result for S — 1, we discuss the S-dependence of the 
phase diagram. The prediction from a bosonization analysis on the decay exponent of the spin 
correlation is verified. 



PACS No: 75.10.Jm, 75.40. Cx 

I. INTRODUCTION 

The study of frustrated quantum spin chains has been 
the subject of much interest for many years. This is 
mainly because these systems exhibit a wide variety of 
exotic phases accompanied with various types of sponta- 
neous symmetry breaking due to the interplay between 
frustration and quantum fluctuations. Among them, one 
of the simplest model may be the quantum spin chain 
with the antiferromagnetic (AF) nearest-neighbor cou- 
pling J\ and the frustrating AF next-nearest-neighbor 
coupling J2. The model Hamiltonian has the form, 
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where Si is the spin-S* operator at site I and A is the ex- 
change anisotropy. Throughout this paper, we consider 
the case of < A < 1 and j = J%/J\ > 0. 

In the classical limit S — > 00, the spin chain (|l|) ex- 
hibits a magnetic long-range order (LRO) in the ground 
state characterized by a certain wavenumber q. The or- 
der parameter is defined by 



m(q) 
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where L is the total number of spins. The LRO is of the 
Neel-type (q = tt) when the frustration is smaller than 
a critical value, j < 1/4, whereas it becomes of helical- 
type for j > 1/4 with a wavenumber q = cos -1 ( — 1 / 4 j ) . 
Both the time-reversal and parity symmetries are bro- 
ken in this helical ordered phase. In the XY-like case 



(0 < A < 1), the helical ordered state possesses a two- 
fold discrete chiral degeneracy characterized by the right - 
and left-handed chirality, in addition to a continuous de- 
generacy associated with the original £7(1) symmetry of 
the XY spin. The chiral order parameter is defined byu 
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Note that this vector chirality ki is distinct from the 
scalar chirality often discussed in the literaturea defined 
by Xi ~ &i-i ■ Si x §1+1: The chirality O k changes its 
sign under the parity operation but is invariant under the 
time-reversal operation, while the scalar chirality changes 
its sign under the both operations. 

In the quantum case S < 00, it seems well established 
that no magnetic LRO emerges at least for < A < 1 in 
contrast to the classical limit. The magnetic LRO (^) is 
destroyed by quantum fluctuations. We note that the ab- 
sence of magnetic LRO in the quantum spin chain (0) is 
proved rigorously for the XY (A = 0) and the Heisenberg 
(A = 1) casesfl By contrast, no theorem prohibiting the 
spontaneous breaking of the parity symmetry has been 
known. Accordingly, there remains a possibility of the 
appearance of a novel "chiral" ordered phase in which 
only the chirality (ra) exhibits a LRO without the mag- 
netic helical LRO (g). This phase breaks only the parity 
symmetry spontaneously with preserving both the time- 
reversal and translational symmetries. 

For the S — 1/2 chain, Nersesyan et al. predicted, 
using the bosonization technique combined with a mean- 
field analysis, that in the XY case (A = 0) the system 
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might exhibit a chiral ordered phase with gapless exci- 
tations for large j.u This gapless chiral phase, however, 
has not been identified in our previous numerical work 
in which the Binder parameter of the chirality was calcu- 
lated numerically for the 5 = 1/2 XY chain with up 
to L =ri20 sites using the exact-diagonalization (ED) 
method.El Although Aligia et alB pointed out that the 
system size L = 20 might be insufficient to deal with 
the chirality in the large j region where the gapless chi- 
ral phase was expected, the question whether the chiral 
phase is realized in the 5 = 1/2 case has not been clari- 
fied so far. 

Meanwhile, the situation seems less controversial in 
the 5 = 1 case. In the previous works, we deter- 
mined the ground-state phase diagram of the 5=1 
XX Z chain by means of both the ED and the density- 
matrix renormalization-group (DMRG) calculations, and 
showed that the gapless chiraL. phase appeared in a 
broad region of the j-A plane.EiL] This observation was 
supported-, by the recent field-theoretical analysis by 
Kolezhukja in which he obtained the gapless chiral phase 
not only for 5=1 but also for general integer 5. Le- 
cheminant et al. extended the bosonization analysis by 
Nersesyan et al. to general 5 and concluded the appear- 
ance of the gapless chiral phase.L 

An interesting observation which was first revealed by 
our numerical study of the 5=1 chain-,is that there 
exist two different types of chiral phasesWi In the one 
of them, the gapless chiral phase, the chiral LRO exists 
while the string and spin correlations decay algebraically 
suggesting gapless excitations. In the other, called the 
gapped chiral phase or the "chiral Haldane" phase, the 
chiral and string LRO's coexist and the spin correlation 
decays exponentially suggesting a finite energy gap. The 
chiral Haldane phase exists in a very narrow but finite re- 
gion between the Haldane and the gapless chiral phases. 
For example, in the pure XY case (A = 0), the sys- 
tem undergoes two successive transitions with increasing 
j, first at j — ~ 0.473 from the Haldane phase to 

the chiral Haldane phase, and then, at j = ~ 0.490 
from the chiral Haldane phase to the gapless chiral phase. 
The next question which arises naturally is whether the 
gapped chiral phase is even realized for general 5^1, 
or it is specific to 5 = 1. In Ref. ||, Kolezhuk suggested 
that the gapped chiral phase existed also for general in- 
teger 5, although the method used was not applicable to 
half-odd integer 5. Lecheminant et al. showed within 
the bosonization analysis that the gapless chiral phase 
realized at larger j might undergo a transition into the 
gapped chiral phase with decreasing jJor any 5 including 
both integer and half-odd integer SB However, both of 
these analyses were based on some approximations, and 
the question whether the gapped chiral phase exists for 
general 5^1 still remains open. 

The aim of this paper is to examine whether the gapless 
and gapped chiral phases exist in the frustrated XXZ 
chain (Q) for general 5. Furthermore, we wish to clarify 



how the chiral ordered phases change their character as 
the spin quantum number 5 increases, to be connected 
to the helical ordered phase realized in the classical limit 
5 — > co. For these purposes, we numerically determine 
the ground-state phase diagrams of the spin chain (Q) in 
the cases of 5 = 1/2, 3/2, and 2. The.method used is the 
same as that in our previous work.Q Using the DMRG 
method, we calculate appropriate correlation functions 
associated with the order parameters characterizing each 
phase, and analyze their long-distance behaviors. By 
combining the obtained results with those of the 5 = 1 
chain determined previously]!! we find that the gapless 
chiral phase appears in the cases of general 5 > 1/2. By 
contrast, the gapped chiral phase has been identified only 
in the cases of integer 5 (5 = 1 and 2): It has not been 
identified within our present numerical accuracy in the 
cases of half-odd integer 5 (5 = 1/2 and 3/2). It also 
turns out that the region of the gapless chiral phase be- 
comes broader as 5 becomes larger, smoothly converging 
to the region of the helical ordered phase in the 5 — > oo 
limit realized at j > 1/4. 

The plan of this paper is as follows. In Sec. II, we ex- 
plain the method used in this paper. Various correlation 
functions associated with each order parameter are intro- 
duced. We show the results of our DMRG calculation in 
Sec. III. The phase diagrams obtained for 5 = 1/2, 3/2, 
and 2 are presented in Sec. Ill A, B, and C, respectively. 
By combining the results with the 5=1 phase diagram 
determined previously,!!] we discuss the 5-dependence of 
the phase diagram in Sec. IV. We also compare our re- 
sult on the decay exponent of the spin correlation with 
the prediction from the bosonization studies .QU Finally, 
our results are summarized in Sec. V. 



II. CORRELATION FUNCTIONS AND 
NUMERICAL METHOD 

In this section, we introduce various correlation func- 
tions and explain our numerical method. 

For half-odd integer 5, it has been known via previous 
studies that there exist twojiistinct phases, the spin fluid 
(SF) and the dimer phases.E3 The SF phase is character- 
ized by gapless excitations above the singlet ground state 
and an algebraic decay of spin correlations. By contrast, 
the dimer phase is characterized by a finite energy gap 
above the doubly degenerate ground states and an expo- 
nential decay of spin correlations. In the dimer phase, 
both the parity and translational symmetries are broken 
spontaneously. The order parameter characterizing the 
dimer phase is given as 
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Tf = (-l)'SfSf +1 . 
For the 5=1/2 and 5 = 3/2 open chains, we calculate 
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the chiral, dimer, and spin correlation functions defined 
by 

C«W = ^r(«io-r/2Kj +r/2>, (5) 

f~IX („\ _ ^ IPX qX ( qX qX 

L/ diinV' ) — ~gi\ D l a -r/2 D U,-r/2+l x \ D l +r /2°l a +r /2+1 

~~ ^+r/2+l^+r/2+2))' (^) 

C s » = ^<Sf _ r/2 5f o+r/2 ) (a = x,z), (7) 

which are associated with the order parameters (||), (|4|), 
and (^]), respectively. The site number l represents the 
center position of the open chains, i.e., Iq = L/2 for even 
r and lo = (L+l)/2 for odd r. We note that the number 
of site L is set to be even. The notation (• • •) represents 
the expectation value at the lowest energy state in the 
subspace of S z total = 0. 

For integer 5, on the other hand, it has been known 
that there generally exist the SF and the Haldane 
phasesJU except for the 5 = lrff&se where only the Hal- 
dane phase appears for A > O.Lrcl The Haldane phase, in 
which the spin correlation decays exponentially with a 
finite energy gap above the singlet ground state, is |Chat- 
acterized by the generalized string order parameteni3'Eij 

Ostr = ^ ( E 4 S k ) S l- ( § ) 

I \k=l ) 

In the S = 1 case, it has also been found that there exists, 
the "double-Haldane" phase in the region of larger j j3Efl 
although its existence has not been clear in the cases of 
5 > 2.113 In order to identify these possible Haldane and 
double-Haldane phases in the 5 = 2 open chain, we cal- 
culate the generalized string correlation function, 

/io+r/2-1 \ 

Cstr(r) = ^2(5f _ r/2 exp I ^ l \ S k 5 f +r/ 2 ): 

\k=la-r/2 J 

(9) 

in addition to the chiral correlation function C K (r) and 
the spin correlation function C"(r) introduced above. 

We calculate those correlation functions for various 
fixed values of A (J) with varying j (A), and estimate the 
transition point j c (A c ) by examining the r-dependence 
of the correlation functions at long distances. We em- 
ploy the. infinite-system DMRG algorithm proposed by 
White. E3 In the calculation of the 5 = 2 chain, in partic- 
ular, the accelerated algorithm proposed by Nishino and 
Okunishi is used.EZI The number of kept states to is up 
to m = 450, m = 260, and m = 260 for the 5 = 1/2, 
3/2, and 2 cases, respectively. Convergence of the data 
with respect to m is checked by consecutively increasing 
to. The truncation error of the DMRG calculation in- 
creases dramatically as j becomes larger. Accordingly, 
in the calculation of larger j , we need to keep more and 



more states in order to achieve the m- convergence of the 
data. Due to this difficulty, our calculation is limited to 
rather small j, i.e., j < 1.6 for S = 1/2 and j < 1 for 
S = 3/2 and 2. In the infinite-system DMRG algorithm, 
the system size is increased by two at each DMRG step 
and the calculation is continued until the i-convergence 
of the data has been attained. We have performed the 
calculation of typically 1000 DMRG steps (corresponding 
to the system with L = 2000 sites) by confirming the L- 
convergence. We can thereby safely avoid the finite-size 
effect arising from the incommensurate character of the 
spin correlation as pointed out by Aligia et alB 



III. NUMERICAL RESULTS 

A. Spin-1/2 case 

We begin with the 5 = 1/2 case. The frustrated XX Z 
spin chain ([!]) for 5 = 1/2 has been studied extensively so 
far. We first review the known properties of the ground- 
state phase diagram. With increasing j, the system un- 
dergoes a Kosterlitz-Thouless (KT) phase transition from 
the SF phase to the dimer phase at j — j^ 2 -'@ , lii The 
critical value j^ 2 ^ has been estimated accurately for 
< A < 1: It runs from j { d 1/2) (A = 0) ~ 0.33 in the 

XY case to j { d 1/2) (A = 1) ~ 0.2411 in the Heisenberg 
case as A increases. E3 Furthermore, the phase diagram 
is divided into two regions according to the nature of 
the short-range spin correlation. The structure factor 
S(q) has a maximum at q = it for j < j^ 2 ^ whereas 
the maximum of S(q) occurs at an incommensurate po- 
sition q < 7r for j > j^ 2 ^ ■ The Lifshitz point j^ 2 ^ 
in the Heisenberg case (A = 1) has been estimated to 

be j^y 2 ^ — 0.52.E3 Meanwhile, the phase diagram in the 
large j region remains largely unclear. In particular, as 
already mentioned in Sec. I, the question whether the 
chiral phase ever exists for larger j still remains contro- 
versial. 

Now, we show our numerical results on the dimer-chiral 
transition. The calculated chiral, dimer, and spin cor- 
relation functions are shown in Fig.[l] (a)-(c) on log-log 
plots for A = and for several typical values of j. As 
can clearly be seen from Figjl] (a), the chiral correlation 

function C K (r) for j > j^ 2 ^ — 1.26 is bent upward at 
larger r suggesting a finite chiral LRO, while C K (r) for 

j < j^ 2 ^ is bent downward suggesting an exponential 
decay of chiral correlations. Although the data around 
the transition point suffer from the truncation error in- 
herent to the DMRG calculation, we can estimate the 
transition point as j^ 2 " 1 = 1.261q'q3 by taking account 
of the m-dependence of the data shown in the figure. 
Thus, we conclude that the chiral ordered phase appears 
even in the 5 = 1/2 chain for j > j^ 2 ^ — 1.26. The re- 
sult is consistent with the prediction of the bosonization 
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it is bent downward for j > j& ■ 
dimer transition point j^J 2 ^ to be j& 



studyfl 

Meanwhile, as shown in Fig|l] (b), the dimer correla- 
tion function C| im (r) for j < jjg — 1.26 is bent up- 
ward for larger r suggesting a finite dimer LRO, whereas 

We estimate the 
1.26 ±0.01. 

We note that, contrary to the chiral correlation function, 
the m-convergence of C% im (r) has almost been attained 
at 77i — 450. Figure [If (c) shows the spin correlation 
function Cf(r) divided by the leading oscillating factor 
cos(Qr), where Q is the wavenumber characterizing the 
incommensurability of CJ (r) in real space. Here, the m- 
convergence of the data has also almost been attained. 
Although the plots are largely scattered, which might be 
attributed to the possible influence of correction terms 
characterized by wavenumbers Q' ^ Q, it is still clearly 
visible in the figure that the spin correlation is bent down- 
ward for j < j^ 2 ^ suggesting an exponential decay, while 
it exhibits a linear behavior for j > j% suggesting a 
power-law decay. From these behaviors of the dimer and 
spin correlations, we conclude that, as j increases, the 
system exhibits a transition from the gapped phase with 
the dimer LRO (j < j^ 2 ^ 1 ) to the gapless phase without 
the dimer LRO {j > j ( c l /2) ). 



The remaining problem is the relation between 



(1/2) 



and j 



,-(1/2) 



Two possibilities seem to be allowed from 



our data, i.e., (i) = or (ii) < j% 2 \ 

If the case (i) is realized, the system undergoes only one 



phase transition at j = j^ 2 ^ 



Jc2^ 2 ^ between the dimer 



phase and the gapless chiral phase with no dimer order. 
If the case (ii) is realized, on the other hand, the sys- 
tem undergoes two successive transitions on increasing j, 
first at j = f rom the dimer phase to the gapped 

chiral phase (or the "chiral dimer" phase) where both 
the dimer and chiral LRO's coexist with gapfull excita- 
tions, and then at j — j^ 2 ^ from the chiral dimer phase 
to the gapless chiral phase. Unfortunately, rather large 
error bars of j^f and prevent us from determin- 

ing which of the cases is realized. However, our result 
suggests that the chiral dimer phase, if it ever exists, 
appears only in a rather narrow region, less than about 
3% of j c ~ 1.26, between the dimer and gapless chiral 
phases. (We note that in the S = 1 XY chain the gapped 
chiral phase exists for 0.473 <j < 0A9, r whose width cor- 
responds to about 3.6% of j c ~ 0.473.Lrl3) Further work 
will be necessary to settle the question whether the chiral 
dimer phase exists in the 5 = 1/2 case. 

Performing the calculations for several fixed values of 
A (J) with varying j (A), we determine the phase bound- 
ary between the dimer and the gapless chiral phases. The 
obtained transition points j^ 2 ^ (A^ 2 "*) are plotted in 
Fig. H together with the SF-dimer transition points j d 
determined by Okamoto and Nomura.E3 As can be seen 
from the phase diagram, the gapless chiral phase does 



exist in the large j region. At least within our present 
numerical precision, on the other hand, we did not find 
any evidence of the gapped chiral (chiral dimer) phase at 
all estimated points between the dimer and the gapless 
chiral phases. 

Finally, we mention the reason why the ED analysis of 
Ref. H on the chain up to L = 20 sites failed to detect the 
chiral ordering in the S — 1/2 XY chain. We now con- 
sider that this failure can be ascribed to the finite-size 
effect discussed in Ref. ^. The wavenumber Q charac- 
terizing the incommensurability decreases rapidly from 
Q = 7r at j = jl to Q = 7r/2 at j — > oo as j increases. 
In the S = 1/2 XY chain, the shift of Q from tt/2 turns 

out to be smaller than 0.037T for j > j^ 2 ' 1 — 1.26. Mean- 
while, in a finite open chain with L sites, the incom- 
mensurability smaller than 2tt/L cannot be taken into 
account because of the condition that the wavefunction 
should vanish on both ends of the chain. As a conse- 
quence, the system size must be larger than L ~ 70 if 
one wishes to detect the shift of Q of order 0.037T. Thus, 
the ED calculation on chains up to L = 20 might fail 
to extract the true asymptotic properties of the chiral 
Binder parameter for j > j^ 2 ^ ■ 



B. Spin-3/2 case 

In this subsection, we present the results for the S = 
3/2 case. The frustrated 5 = 3/2 chap. (|J) was stud- 
ied for the Heisenberg case (A = 1).BB It was shown 

there that the system was in the SF phase for j < j d 

(3/21 

while it was in the dimer phase for j > j\' . The phase 



(3/2) 



0.33. 



transition of the KT type occurs at j = j d 

(3/2) 

The Lifshita-point was estimated to be at j = j L ( A = 
1) = 0.388.B Besides, the appearance of the gapped and 
gapless chiral phases was suggested recently for the XY 
case (A = 0) by the bosonization study.0 

In Fig. ||, we show the DMRG data of the spin cor- 
relation function C^{r) for A = 0.6 and for several typ- 
ical values of j. As can be seen from the figure, the 
spin correlation changes its behavior as j increases from 

(3/2) 

a power-law decay (j < j d ~ 0.335) to an exponential 

( 3 /2) 

decay (j > j d )• We estimate the transition point be- 
tween the gapless SF phase and the gapped dimer phase 
as jf /2) (A = 0.6) = 0.335 ± 0.015. By interpolating the 
transition points estimated in this way for various A, we 
determine the SF-dimer transition line. The line runs 
from j ( d /2) (A = 0) = 0.334 ± 0.004 in the XY case to 

(3 /2) 

j d (A = 1) = 0.335 ±0.015 in the Heisenberg case, the 
latter being consistent with the estimate by Ziman and 
Schulz, jf /2] (A = 1)~ 0.330 

Let us next consider the dimer-chiral transition. Fig- 
ure H (a) shows the chiral correlation function C K (r) for 
A = 0.6 and for several typical values of j. As shown in 



the figure, C K {r) is bent upward for j > j cl 



(3/2) 



0.404 
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suggesting a finite LRO, while it is bent downward for 

(3 /2) 

j < Jci suggesting an exponential decay. Hence, the 
chiral ordered phase is realized also in the 5 = 3/2 case. 
The transition point at which the chiral LRO sets in is es- 
timated to be j { c l /2) = 0.404±g;{$§. The spin correlation 
function Cjf (r) for A = 0.6 is shown in Fig. [| (b) where 
the data are divided by the oscillating factor cos(Qr). As 
in the S = 1/2 case, the plots are largely scattered sug- 
gesting the existence of non-negligible correction terms. 
Nevertheless, it can be seen that the behavior of Cf(r) 
changes as j increases from an exponential decay to an 

(3 /2) 

algebraic decay at around j = j c2 = 0.410 ±0.010. We 

(3 /2) 

thus conclude that the system for j > ' is in the gap- 
less chiral phase. Unfortunately the error of our estimate 
of j^ 2 ^ is quite large. This large error is ascribed to the 
poor m-convergence of the data of Cf(r) shown in Fig. 
[| (b) and of the data of the dimer correlation function 
C| im (r) to be mentioned below. (Note that the situation 
here is different from the 5 = 1/2 case where the m- 
convergence of the spin and dimer correlation functions 
is almost achieved at m = 450: See Fig. |l].) Because of 
the errors in the estimates of and j% , we can not 

(3/2) . . " 

determine whether j c2 is either larger than or equal to 
^(3/2) rp nug ^ j n Qur p resen ^ calculation, the gapped chi- 
ral phase has not been identified in the 5 = 3/2 case 
as well as in the 5 = 1/2 case. We have performed the 
same calculation for various A with varying j, but have 
not been able to identify the gapped chiral phase at any 
j and A. 

In Fig. |^ (a), we show the dimer correlation function 
C^ im (r) for A = 0.6 and for several typical values of j. At 
a glance, the data may look like indicating a finite dimer 
LRO even in the SF phase (j < jf /2) ~ 0.335) and in 

(3/2) 

the gapless chiral phase (j > j c2 — 0.410) where no 
dimer LRO is to be expected. The dimer correlations 
in these phases, however, rapidly decrease with increas- 
ing m, whereas the m-convergence has almost been at- 
tained in the intermediate range of j corresponding to 

(3/2) (3/2) 

the gapped phase (j d < j < j c2 )• To elucidate the 
m-dependence of the long-distance value of the dimer 
correlation, we plot in Fig. [5] (b) the value of C^ im (r = 
100), which is .expected to be a good approximation of 
C£ m (r -> oo).E3 The figure shows that C| im (r = 100) 
in the gapped phase converges to a nonzero value as 
vn — ► 00, while it decreases toward zero in the SF and 

the gapless chiral phases. Thus, we infer that the appar- 

(3/2) 

ent finite dimer LRO observed in Cj im (r) for j < j d 

(3 /2) 

and jj-> y c2 is a spurious effect due to the truncation 
errorEa The variation of the data with varying j near 
the expected phase boundaries is rather gentle, which 
prevents us from estimating the transition points accu- 
rately from Fig. || (b). Nevertheless, we may conclude 

that the system in the gapped region j^ 2 * 1 < j < jit 
possesses a true dimer LRO and is indeed in the dimer 
phase. 



We show in Fig. |6| the obtained ground-state phase di- 
agram of the 5 = 3/2 chain including the SF-dimer and 

the dimer-chiral transition lines, and j^i '. We 

note that the dimer phase exists even in the XY case 

(A = 0) for j ( d 3/2} ~ 0.334 < j < j^ /2) ~ 0.339 although 

(3/2) 

its width is quite narrow. The value of j cl becomes 
larger as A becomes larger, and the dimer-chiral phase 
boundary appears to tend to the Heisenberg line A = 1 
as j — > 00. 

C. Spin-2 case 

In this subsection, we present our results for the 5 = 2 
case. It has been known that the 5 = 2 chain with only 
the nearest-neighbor coupling exhibits a phase transition 
at A = A^ (j =. 0) ~ 0.966 between the SF and the 

Haldane phascs.EJ Meanwhile, the Lifshitz point in the 

(2) 

Heisenberg case was estimated to be at j — j L (A = 
1) = 0.325.E£l Very recently, the existence of the gap- 
less and gapped chiral phases has been suggested-by the 
large-5 approacbtl and the bosonization method.u 

Let us first consider the SF-Haldane transition at 
j = jjj . In Fig. we show the generalized string and 
spin correlation functions C s t T (r) and Cf (r) for A = 
and for several typical values of j on log-log plots. It 

can be seen in Fig. R (a) that C s t r (V) is bent upward 

(2) 

for j > jir 1 ~ 0.280 suggesting a finite LRO, while it 

(2) 

shows a linear behavior for j < j H suggesting an alge- 
braic decay. Meanwhile, as shown in Fig. fj] (b), (r) 

(2) 

decays exponentially for j > j H indicating a finite gap, 

(2) 

whereas it decays algebraically for j < j H '. We thus 

(2) 

conclude that there occurs a phase transition at j = j H 
between the SF and the Haldane phases. The critical 
value is estimated to be = 0.280±0.006. We esti- 
mate the critical points for various fixed A (J), and deter- 
mine the phase boundary. The transition line smoothly 
connects the point j = J^(A = 0) = 0.280 ± 0.006 in 
the XY case to the point of the case of no frustration, 
A = A (2) (j = 0) = 0.96 ± 0.01. (See the phase dia- 
gram shown in Fig. |[) The latter estimate is consistent 
with the previous estimate by Nomura and Kitazawa, 
A {2 \j =0) ~ 0.966.0 

Next, we consider the transition between the Haldane 
and the chiral phases. Fig. || (a) exhibits the chiral cor- 
relation function C K (r) for A = and for several typical 

values of j. As can clearly be seen from the figure, C K (r) 

(2) 

exhibits a finite LRO for j > j cl ~ 0.318 whereas it 

(2) 

exhibits an exponential decay for j < j cl . We estimate 
the transition point where the chiral LRO sets in to be 
jif = 0.318 ± 0.001. As shown in Fig. | (b), the gen- 
eralized string correlation C s t r (0 exhibits a finite LRO 

(2) 

for j < y c2 ~ 0.324 whereas it decays algebraically for 
j > . Meanwhile, as shown in Fig. || (c), the spin 
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correlation function Cjf (r) divided by the leading oscil- 

(21 

lating factor cos(Qr) decays exponentially for j < j c2 

(2) 

whereas it decays algebraically for j > j c2 . From these 
observations, we estimate the transition point where the 
excitation spectrum becomes gapless and the generalized 
string order vanishes to be jf 2 = 0-324to'go2- Here the 

(2) 

estimate of j c2 is quite close to, but distinctly larger than 

that of jj? = 0.318±0.001. Indeed, the data of j = 0.320 
in Fig. @ clearly show the existence of an intermediate 
phase, the gapped chiral phase, where the chiral and the 
generalized string LRO's coexist and the spin correla- 
tion decays exponentially. Hence, we conclude that, as 
in the case of the 5 = 1 chain, the gapped chiral (chiral 
Haldane) phase exists also in the 5 = 2 chain in a very 

(21 (21 

narrow but finite region (j^ < j < j&j between the 

(2l 

Haldane phase {j < y cl r ) and the gapless chiral phase 

(2) 

(j > j c2 ). Performing the same calculations for various 

(2) (2) 

fixed A, we estimate the transition points j cl and j c2 . 
The chiral Haldane phase is also found in a certain range 
of j for A = 0.2, 0.4, 0.6, and 0.8. By contrast, we can 
not confirm its existence for A > 0.8 corresponding to 
j > 0.5 within our numerical accuracy. This is due to the 
large truncation error growing drastically with increas- 
ing j, which prevents us from precisely determining the 
transition points. 

The obtained phase diagram of the 5 = 2 chain is 
shown in Fig. ^. The gapless chiral phase appears in a 
quite broad region of large j. The boundary of the re- 

(2) 

gion with the chiral LRO, j cl , rises steeply as j increases 
from the point on the XY line, j^f (A = 0) ~ 0.318, 
tending to the Heisenberg line A = 1. We note that, 
although our estimates of the points where the chiral 
LRO sets in are very close to A = 1 for j > 0.6, i.e., 
A^ } = 0.9975 ± 0.0025 for all j = 0.6, 0.7, 0.8, and 1.0, 
the chiral LRO is not observed on the Heisenberg line 
A = 1 where there no longer exists the two-fold discrete 
chiral degeneracy. 

We finally refer to our numerical results on the possi- 
ble double-Haldane (DH) phase. The DH phase was first 
found in the-large j region of the frustrated 5 = 1 Heisen- 
berg chain. til In the DH phase, the next-nearest-neighbor 
coupling Ji is dominant and the system can be regarded 
as two Haldane subchains coupled by the weak inter- 
subchain coupling J\. The DH phase is characterized 
by the absence of the string, spin, and chiral LRO's.E3 It 
has been shown that in the frustrated 5=1 XX Z chain 
with A > 0.95, there occurs a first order phases-transi- 
tion between the Haldane and the DH phases. Lili-a The 
string order parameter vanishes discontinuously at the 
transition point. In our present study on the frustrated 
5 = 2 Heisenberg chain, such a vanishing of the string 
order has not been observed: The generalized string cor- 
relation function C s t r (r) in the Heisenberg case (A = 1) 
exhibits a finite LRO for an entire region studied here, 
< j < 1, while both the spin and chiral correlation 



functions CJ (r) and C K (r) show an exponential decay 
for < j < 1 suggesting a gapfull excitation. We note 
that, as j increases, the extrapolated value C s t r (V — » oo) 
decreases rapidly around j ~ 0.4. Although the rapid 
drop of the string LRO might be a sign of a phase transi- 
tion from the Haldane phase to a new intermediate phase 
between the Haldane and the DH phases as suggested in 
Ref. [l5|, the situation here remains largely unclear. Fur- 
ther work is required to clarify the details of the transi- 
tion at and near A = 1 and the possible existence of the 
DH phase for j > 1. 



IV. S-DEPENDENCE 

Based on our present results foi^.5 = 1/2, 3/2, and 2 
and our previous result for 5 = 1,0 we now discuss how 
the ground-state phase diagrams change as 5 increases 
from 5=1/2 toward the classical limit 5 — > oo. We are 
interested particularly in the way how the classical limit 
is achieved from the quantum phases analyzed above. 

By comparing the obtained phase diagrams for 5 = 
1/2, 1, 3/2, and 2, we deduce several features of the phase 
diagrams. The first is about the fate of the dimer and the 
Haldane phases. These phases are pure quantum ones in 
the sense that singlet spin pairs play an essential role in 
stabilizing them. Hence, it is natural to expect that the 
regions of these quantum phases become narrower as 5 
gets larger. This feature can clearly be seen in the phase 
diagrams obtained for 1/2 < 5 < 2: The dimer phase 
for 5 = 3/2 is narrower than that for 5 = 1/2, while the 
Haldane phase for 5 = 2 is narrower than that for 5 = 1 . 
It thus seems reasonable to expect that the dimer and the 
Haldane phases continue to become narrower for 5 > 2, 
and eventually vanish in the classical limit 5 — > oo. 

The second feature concerns the chiral ordered phases. 
Recently, the 5-dependence of the boundary of the chi- 
ral ordered phases,-was examined via the field-theoretical 
large-5 approach.B It was shown there that, as 5 in- 
creased, the region of the gapless chiral phase converged 
smoothly to that of the helical ordered phase in the clas- 
sical limit 5 — > oo, whereas the gapped chiral phase van- 
ished asymptotically. This feature of the chiral phases 
might be understood intuitively by considering the role of 
quantum fluctuations in the symmetry breaking. In the 
helical ordered state with a finite magnetic LRO realized 
in the classical limit 5 — > oo, both the discrete Zi par- 
ity and the continuous f/(l) spin symmetries are broken 
spontaneously. In the gapless chiral phase, quantum fluc- 
tuations marginally recover the continuous U(l) symme- 
try yielding the quasi-long-range spin order, with keeping 
the discrete Zi parity symmetry broken. In this sense, 
the gapless chiral phases can be regarded as a quantum 
remnant of the classical helical phase. By contrast, the 
gapped chiral phase should be regarded as a pure quan- 
tum phase since it exhibits a topological LRO, i.e., the 
string LRO, in the chiral Haldane phase for integer 5. (In 
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the chiral dimer phase for half-odd integer 5, if any, the 
topological LRO is the dimer LRO, although we donot 
find evidence of such a phase in the present study.) We 
therefore expect that the gapped chiral phase should van- 
ish in the classical limit S — ► oo. The observed behaviors 
of the chiral phases in the obtained phase diagrams are 
consistent with the above expectation. With increasing 
5, the chiral transition line j'^f (A) approaches that of 
the helical ordered phase, j = 1/4. Meanwhile, the region 
of the chiral Haldane phase found for integer 5 shrinks 
as 5 increases. As an example, we list in Table | the 

( S) ( S) 

estimated values of j cl and j c2 m the XY case (A = 0) 

and their ratio (j^ — Jci^VJc? f° r integer 5, which is 
a measure of the relative stability of the chiral Haldane 
phase. 

The third feature concerns the SF phase. One may nat- 
urally expect that, as 5 increases, the region of the SF 
phase converges to that of the Neel phase in the classical 
limit 5 — > oo. Somewhat unexpectedly, in the obtained 
phase diagrams, the SF phase grows as 5 becomes larger, 
exceeding the classical phase boundary of the Neel phase, 
j = 1/4. (See the phase diagrams and Table | in which 
the estimates of j d and j H in the XY case are listed.) 
We consider that this "overshooting" is due to the rapid 
shrink of the dimer and the Haldane phases and that the 
boundary of the SF phase eventually "turns back" con- 
verging to the classical phase boundary j = 1/4 for large 
enough 5. In order to confirm this conjecture, the anal- 
ysis of larger 5 is needed, which is beyond the scope of 
the present work. 

Finally, we wish to compare our numerical result on 
the decay exponent of the spin correlation with the pre- 
diction from the bosonization analysis Applying the 
bosonization technique and the mean-field approxima- 
tion, the authors of Refs. |^ and ^ predicted that for large 
enough j the spin correlation function Cf (r) in the XY 
case should exhibit the asymptotic behavior, 

C%{r)=Acos(Qr)r- r '*, (10) 

where A is a numerical constant and the exponent r\ x is 
given by r\ x — 1/(85). We fit our numerical data for the 
5 = 1/2, 3/2, and 2 XY chains with Eq. @, taking A, 
Q, and rj x as fitting parameters. Resulting estimates of 
the decay exponent rj x are shown in Fig. |l0| together with 
our previous estimate for the 5=1 XY chain.0 As can be 
seen from Fig. [r], the estimates of r\ x decrease monoton- 
ically as j increases. For 5 > 1, the asymptotic j ' —> oo 
value has almost been attained at around j — 1, which is 
in good agreement with the predicted value 1/(85). For 
the 5 = 1/2 case, on the other hand, we can not reach 
the asymptotic j ' —>■ oo regime: The estimated r] x still 
continues to decrease even around j ~ 1.5, the largest j 
value for which we can get reliable data. The estimated 
r) x , however, shows a tendency to further decrease toward 
the predicted value r] x — 1/4 as j increases, which is con- 
sistent with the result of the bosonization study. Thus, 



our result can be considered as a numerical support of 
the bosonization analysis. 



V. SUMMARY 

In the present work, we have studied the ground-state 
properties of the frustrated spin-5 XX Z chains (|l|), es- 
pecially paying attention to the chiral ordered phases in 
which only the chirality exhibits a finite LRO without 
the standard magnetic LRO. We have used the infinite- 
system DMRG method to calculate the correlation func- 
tions associated with the spin, chiral, dimer, and string 
order parameters. By analyzing the long-distance behav- 
ior of the correlation functions, we have determined the 
ground-state phase diagrams of the 5 = 1/2, 3/2, and 2 
chains for < A < 1 and j > (Figs, g § and |). 

By comparing the obtained results with our previous 
result for the 5 = 1 chain, we reach the following pic- 
ture of the ground-state phase diagram. In the integer 5 
chains (5 = 1 and 2), the phase diagram consists of four 
different phases, i.e., the SF, Haldane, gapless chiral, and 
gapped chiral (chiral Haldane) phases. In the half-odd in- 
teger 5 chains (5 = 1/2 and 3/2), on the other hand, we 
have found three phases, i.e., the SF, dimer, and gapless 
chiral phases. It thus turns out that the gapless chi- 
ral phase appears for general 5 > 1/2. For integer 5, 
the gapped chiral (chiral Haldane) phase exists in a nar- 
row region between the Haldane and the gapless chiral 
phases. For half-odd integer 5, by contrast, the gapped 
chiral (chiral dimer) phase has not been identified within 
our numerical precision: The rather large truncation er- 
ror of the DMRG calculation prevents us from verifying 
whether the chiral dimer phase exists or not. Our results 
suggests, however, that the chiral dimer phase, if it ever 
exists, appears in a narrow region between the dimer and 
the gapless chiral phase. Further work will be required 
to solve the problem. 

We have also discussed the 5-dependence of the phase 
diagrams. The obtained phase diagrams indicate that, 
as 5 increases toward the classical limit 5 — > oo, the 
region of the gapless chiral phase converge smoothly to- 
ward that of the classical helical phase, while the pure 
quantum phases, i.e., the dimer, Haldane, and gapped 
chiral phases, become narrower and eventually vanish in 
the 5 — > oo limit. The prediction from the bosonization 
study that the decay exponent of the spin correlation in 
the XY case (A = 0) takes a value rj x — 1/(85) in the 
j ' —>■ oo limit has been verified. 

Finally, we wish to touch upon the possible experi- 
mental realization of the gapless chiral phase in a quasi- 
one-dimensional compound CaV2 04. This material is 
expected to be described by the frustrated 5 = 1 chain 
(Q) where the AF next-nearest-neighbor coupling J2 is 
comparable to the AF nearest-neighbor coupling Ji . Re- 
cently, Kikuchi made measurements on magnetic suscep- 
tibility and on 51 V NMR, and showed that the system 



7 



had gapless excitations above the ground stateB For the 
frustrated S — 1 chain (Q), theoretical studies indicate 
that there is no gapless phase for j > and A > ex- 
cept for the gapless chiral phase. Hence, CaV204 might 
be a promising candidate for the realization of the chiral 
ordered phase. For the future, it might be interesting 
to calculate thermodynamic properties at finite temper- 
atures to compare them with the experimental data on 
this compound. 
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TABLE I. The estimates of jf, j^', and j£ in the XY case (A = 0) for half-odd integer S (S = 1/2 and 3/2) and 
for integer S (S = 1 and 2). The ratio (j^f — icl')/iii' f° r integer S 1 is also listed. 



5 


if 




AS) 

Jcl 


AS) 




1/2 


0.33 a 




1 26 +u ' ul 

1 - AO -0.03 


1.26 ±0.01 




3/2 


0.334 ± 0.004 




0.339 ± 0.001 


UJ ™-0 .002 




1 







0.473 ± 0.001 


u.y:»u_ 005 


0.036 


2 




0.280 ± 0.006 


0.318 ±0.001 


n o 24 +0.006 
u -°^-o .002 


0.019 



a Rcf. |19|. 

The data for 5 = 1 are from Ref. 0. 
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(c) spin correlation S=1/2, A=0 




FIG. 1 . Various correlation functions of the S =1/2 chain for A = and for several typical values of j around the dimer-chiral 
transition: (a) chiral correlation function C K (r); (b) dimer correlation function ( — l) r Cj im (r); (c) spin correlation function CJ (r) 
divided by the oscillating factor cos(Qr). The number of kept states is m = 450. To illustrate the m-dependence, we also 
indicate by crosses the data with m = 400 and 350 for j — 1.24, 1.26, and 1.28 in figure (a) where the m-dependence is relatively 
large. In other cases, the truncation errors are smaller than the symbols. 
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FIG. 2. The ground-state phase diagram of the S = 1/2 chain, where j and A denote the ratio J2/J1 and the exchange 
anisotropy, respectively, defined in Eq. (|l|). The diamonds and circles represent the transition points j^ 2 ^ and j^/ 2 ^ ■ 



S=3/2,A=0.6 




FIG. 3. The spin correlation function ( — l) r Cf (r) of the S = 3/2 chain for A = 0.6 and for several typical values of j around 
the SF-dimer transition. The number of kept states is m — 220 for all j. The truncation errors are smaller than the symbols. 
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FIG. 4. Various correlation functions of the S = 3/2 chain for A = 0.6 and for several typical values of j around the 
dimer-chiral transition: (a) chiral correlation function C K (r); (b) spin correlation function C% (r) divided by the oscillating 
factor cos(Qr). The number of kept states is m = 260. We also indicate by crosses the data with m = 220 and 180 for j = 0.400 
and 0.404 in figure (a), and the data with m = 220 for j = 0.400,0.404,0.406, and 0.410 in figure (b). In other cases, the 
truncation errors are smaller than the symbols. 



(a) 



S=3/2, A=0.6 



[] □ i 



„ „ n nnn rTTi i 1 1 1 1 1 mi i i i i 1 1 1 1 1 iiim i i i 1 1 1 1 iittI 



10" 



o 



10" 



x xxxxxxx»o«s«xxxx>ooooo« 



X X 

ii x x xx 

^^ ++ XXxx xxxXKXxxxxxxxxx,>»«i<xxxxxxx«o» 



^ . i -~vwv* x X XXXXXXXXXWX XXXXXXXXW 

^M^vfettt 1- i 



• : j = 0.320 ■ : j = 0.360 a : j = 0.400. 
o : j = 0.340 □ : j = 0.380 a : j = 0.420 



S=3/2, A=0.6 



10' 



10' 




FIG. 5. (a) The r-dependence of the dimer correlation function C| im (r) of the S = 3/2 chain for A = 0.6 and for several 
typical values of j. The number of kept states is m = 220. To illustrate the m-dependence of the data, we also plot the data 
with m = 180 and 140 by crosses for j — 0.320, 0.340, and 0.420. The truncation errors for the other j's are smaller than the 
symbols, (b) The j-dependence of the value of Cj im (r = 100) for various m. The dotted lines represent the phase boundaries, 



j = jf /2) (A = 0.6) ~ 0.335 and j = j { J, /z '(A = 0.6) ~ 0.410 



(3/2), 
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S=3/2 




FIG. 6. The ground-state phase diagram of the S = 3/2 chain, where j and A denote the ratio J2/J1 and the exchange 



anisotropy, respectively, defined in Eq. (111). The diamonds and circles represent the transition points j^ 2 ^ and j£ 
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(b) spin correlation 
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FIG. 7. Various correlation functions of the S — 2 chain for 
transition: (a) generalized string correlation function C s tr{r); 
states is m — 220. The truncation errors are smaller than the 



A = and for several typical values of j around the SF-Haldane 
(b) spin correlation function ( — l) r Cs (r). The number of kept 
symbols. 
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(a) chiral correlation S=2,A=0 
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(b) string correlation S=2,A=0 
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(c) spin correlation S=2,A=0 
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FIG. 8. Various correlation functions of the S — 2 chain for A = and for several typical values of j around the Haldane-chiral 
transition: (a) chiral correlation function C K (r); (b) generalized string correlation function C s tr(r); (c) spin correlation function 
Cs(r) divided by the oscillating factor cos(Qr). The number of kept states is m = 260. To illustrate the m-dependence, we also 
indicate by crosses the data with m = 220 and 180 for j — 0.318 in figure (a), j — 0.324 and 0.330 in figure (b), and j — 0.324 
in figure (c). In other cases, the truncation errors are smaller than the symbols. 
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S=2 




FIG. 9. The ground-state phase diagram of the S = 2 chain, where j and A denote the ratio J2/J1 and the exchange 
anisotropy, respectively, defined in Eq. (|l|). The diamonds, circles, and squares represent the transition points , and 

, respectively. 
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FIG. 10. The estimated decay exponent rj x as a fiuaction of j for the S = 1/2, 1, 3/2, and 2 XY chains. The dotted lines 
represent the prediction of the bosonization analysisu'u at j — > 00, T] x = 1/(85'). 
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